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Abstract 



The high density interior of a neutron star is expected to contain supercon- 
ducting protons and superfluid neutrons. Theoretical estimates suggest that the 
protons will form a type II superconductor in which the stellar magnetic field 
is carried by flux tubes. The strong interaction between the flux tubes and the 
neutron rotational vortices could lead to strong 'pinning', i.e. vortex motion 
could be impeded. This has important implications especially for pulsar glitch 
models as it would lead to a large part of the vorticity of the star being decoupled 
from the 'normal' component, to which the electromagnetic emission is locked. 
In this paper we explore the consequences of strong pinning in the core on the 
'snowplow' model for pulsar glitches (Pizzochero 2011), making use of realistic 
equations of state and relativistic background models for the neutron star. We 
find that in general a large fraction of pinned vorticity in the core is not com- 
patible with observations of giant glitches in the Vela pulsar. The conclusion is 
thus that either most of the core is in a type I superconducting state or that the 
interaction between vortices and flux tubes is weaker than previously assumed. 
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Introduction 



Neutron Stars (NSs) allow us to probe the state of matter in some of the most extreme 
conditions in the universe. Not only can the density in the interior of these very compact 
objects exceed nuclear saturation density, but NSs also host some of the strongest magnetic 
fields in nature, with intensities of up to ~ 10 15 G for magnetars. Not surprisingly modelling 
such complex objects requires the use of some poorly understood physics. 

In particular the star will rapidly cool below the critical temperature for the neutrons 
to be superfluid and the protons to be supe rconducting. The protons of the outer core are 
predicted to form a type II superconductor (jMigdallll959l ; iBaym. Pethick fc Pined 1 19691 ). in 
which the magnetic flux is confined to flux tubes, inside which the magnetic field strength 
is of the order of the lower critical field for superconductivity, B c w 10 15 G. However above 
a critical density of approximately p c ~ 3 x 10 14 g/cm 3 one expects a transition to type I 
superconductivity, in which the formation of fluxtubes is no longer fayourab le but rather 
the magnetic field is contained in regions of normal protons ( jSedrakianl 120051 ) . Given that 
the critical density for this transition is easily reached in NS interiors it is possible that a 
sizeable portion of the star may in fact be in a type I superconducting state ( Uonesl 120061 ). 



The dynamics of the outer core plays a crucial role in the interpretation of various astro- 
physical phenomena, such as pulsar glitches, timing noise, precession and fluid oscillations. 
In particular pulsar glitches are sudden increases in the otherwise steadily decreasing rota- 
tional frequency of a pulsar. Although their origin is still debated it is generally thought 
that these phenomena are the direct manifestation of a superfluid component inside the star, 
which is only weakly coupled to the normal component due to the interaction between the 
quantized neutron vortex lines and the charged particles in the crust or core. In particular if 
vortices can "pin" (i.e. are strongly attracted) to the Coulomb lattice in the crust, they can 
decouple the neutron superfluid from the normal component (to which the electromagnetic 
emission is locked) and t heir sudden depinning will give rise to a rapid transfer of angular 
momentum, i.e. a glitch (jAnderson fc Itoh 1975 ; Alpar 1977 ; Pines et al. 1980 ; Alpar et al. 



1 98 ll : I Anderson et al.l 119821 ) . Recent work has shown that thi s scenario can successfully ac 



count for the distribution in glitch sizes and waiting times (jWarszawski fc Melatod 12008 



Melatos &; Warszawskil 120091 ; IWarszawski &; Melatosl 1201 1. 2012) and describe the size and 



relaxa tion timescales of giant glitches in the Vela pulsar ( IPizzocherol 1201 ll ; lHaskell et al. 
2012ah - 



An important issue to address is, however, whether vortices will only pin to the crustal 
lattice or whe ther they a re pinned to flux tubes if the outer core is in a type II supercon- 
ducting state (jLinki 120031 ) . thus effectively decoupling a large fraction of the stellar moment 
of inertia from the crust. Furthermore if vortices are pinned in the core this is likely to lead 
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to the onset of turbulence and may play an important role in pulsar 'timing noise' ( ILink 
2012bl ). The interaction between flux tubes and yortices can also have a strong impa ct on 
the gravitational wave driven r-mode instability ( IHo et al.ll201ll ; lHaskell et al.ll2012bf ) and 
on NS precession (ILinkl 120031 ). 



In this Let ter we inves t igate the effect of vortex pinning in the core on the "snowplow" 
glitch model of iPizzocherol (120111). We extend the model to realistic equations of state and 
relativistic stellar models, as in ISeveso et al~ ( 2012a ) and show that, in general, one cannot fit 
the size and postglitch jumps in frequency derivative of Vela giant glitches if a large portion 
of the core vortices are pinned. This points towards the fact that most of the core could in 
fact be in a type I superconducting state, or that the vortex/flux t ube interaction is weaker 
than previously assumed, as some microphysical estimates suggest (lBabaevll2009l ). 



Furthermore iGlampedakis fc Anderssonl (120111 ) recently showed that vortex pinning in 
the core is likely to be a short lived phenomenon that may only be relevant in a short 
period of a NS's life and in m agnet ars, and in their hydrodynamical model of giant pulsar 
glitches lHaskell et al.l (12012al ) also find that vortex pinning in the core is inconsistent with 
the observed post-glitch relaxation timescales in the Vela. 



2. The "snowplow" model 



The starting point of our investigation will be the "snowplow" model for glitches of 
Pizzocherol (120111 ). which we briefly review here. We take the NS to be a two component 
system, where one of the components, the so-called "normal" component, is given by the 
crust and all charged components tightly coupled to it by the magnetic field. The other, 
the "superfluid" , is given by the superfluid neutrons in the core and crust. The superfluid 
rotates by forming an array of quantized vortices which carry the circulation and mediate an 
interaction between the two components know n as Mutual Friction, which in th e core can 
couple the two fluids on timescales of seconds (lAndersson. Sidery & Corner! 120061). Vortices 



can, however, also be pinned to ions in the crust or f l ux-tubes in the core (lAnderson fc Itoh 



1975 




1998 





Link 



Alpar 119771 ; iPines et al.lll980t lAlpar et al.lll98ll ; lAnderson et al.lll982t iRuderman et al. 



20031 ) . As a consequence vortex motion is impeded and the superfluid cannot 
spin-down, effectively decoupling it from the normal component which is spinning down 
due to electromagnetic emission. If a lag builds up between the superfluid and the normal 
component this will, however, give rise to a Magnus force acting on the vortices, which takes 
the form f m = K>p s Cl x (v„ — v s ), where f m is the force per unit length, k = 1.99 x 10~ 3 cm 2 s _1 
is the quantum of circulation, p s is the superfluid density, O is the unit vector pointing along 
the rotation axis, v„ is the velocity of the vortex lines and v s is the velocity of the superfluid. 
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We assume that the neutrons are supe rfluid thro u ghout the star and take p s = (l — x p )p, with 



J r. 



the proton fraction calculated by IZuo et al.l (120041 ) . Once the Magnus force integrated 



over a vortex exceeds the pinning force, the vortex will unpin and be free to move out. 



We follow the procedure of ISeveso et al.l (j2012al ) and int egrate the relativistic equa tions 
of stel lar structure for two realistic equat ions of state, SLy (IDouchin fc Haensell 120011 ) and 
GM1 (IGlendenning &: Moszokowskil 119911 ). We assume straight vortices that cross through 
the core (IZhou et al.ll2004h an d for the pinning for c e per unit length f p we use the realistic 
estimates of iGrill et al.l ( 120121 ); iGrill fc Pizzocherol (120121 ) . Balancing the pinning force to 
the Magnus force and integrating over the vortex length allows us to calculate the lag at 
which the vortices will unpin in different regions. The normalization of f p is chosen in such 
a way as to give an inter-glitch waiting time T g = Af2 max /|f2| of approximately 2.8 yrs for 
the Vela pulsar, where Af2 max is the maximum of the critical unpinning lag. 

If there is no pinning in the core vortices will unpin and move out toward the crust, where 
they encounter a steeply increasing pinning potential and repin. This leads to the creation of 
a thin vortex sheet that moves towards the peak of the potential, the so-called "snowplow" 
effect. Once the maximum of the critical lag has been reached the vortices can no longer 
be held in place and the excess vorticity is released catastrophical ly, exchanging an gular 
momentum with the normal component and giving rise to a glitch ( |Pizzocheroll201ll ). We 
assume that this is the mechanism that gives rise to giant glitches, i.e. glitch es with steps in 
the sp in rate AQ g i « 10 -4 rad that are observed in the Vela and other pulsars ( lEspinoza et al. 



20111). Smaller glitches are lik e ly to be triggered by crust quakes or random vortex avalanc hes 
(TWarszawski fc Melatoslbood : iMelatos fc Warszawskilbood : [warszawski fe Melatosl boill ) . 



We can easily calculate the number of vortices in the vortex sheet once it has reached the 
peak of the potential as N v = ^r^ ax Af2 max , where r max is the cylindrical radius at which the 
maximum of the critical lag is located and Af2 max the value of said maximum. The angular 
momentum exchanged as the vortices move out and annihilate is then given by: 



AL a 



2kQN v 



xdx 



l(x)/2 



p(V x 2 + z 2 )dz 



(1) 



where Q = I s / hot is the superfluid fraction of the moment of inertia, R c is the radius of the 
inner crust where the vortices annihilate (taken at neutron drip density), l{x) is the length of 
a vortex at a given cylindrical radius x and p is the density. The glitch observables can then 
be derived as AQ gl = AL// tot [l - Q(l - Y gl )} and Af^oo = [Q(l - Y gl )]/{1 - Q(l - Y gl )], 
where AQ g i is the step in angular velocity due to the glitch and AQgi/Qoo is the instantaneous 
step in the spin-down rate immediately after the glitch, relative to the steady state pre-glitch 
spindown rate Cl^. We have also introduced the parameter Y g \ which represents the fraction 
of superfluid moment of inertia which is coupled to the crust during the glitch. Given that 
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the rise time r r is very short (less than a minute (IDodson. McCulloch fc Lewis! |2002| )) it 
is likely that only a small fraction of the core will be coupled to the crust on this short 
timescale, with the rest of the star recoupling gradually on longer timescales and giving rise 
to the observed exponential post-glitch relaxation (see lHaskell et al.l (j2012al ) for a detailed 
di scussion of this issue). The best ob servational upper limits on the rise time are r r < 40 
s (IDodson. McCulloch fc Lewis! 120021 ) from the Vela 2000 glitch, while an interesting lower 
limit of t > 10~ 4 ms can be derived from the non detection of a GW signal from the Vela 
2006 glitch (IWarszawski fc Melatoa 120121 ). Theoretical estimates, on the other hand, give 
r r w 1 — 10 s (lHaskell et al.ll2012al ). which thus easily allow for the angular momentum in 
(JTJ) to be exchanged during the short rise times observed in radio pulsars. 

Let us now consider the motion of a vortex if the NS core is a type II superconductor. 
As a vortex approaches a flux tube its magnetic energy will increase if they are alligned or 
decreas e if they are antialigned reasulting in an energy per intersection of approximately E p « 
5 MeV ( iRuderman et al.lll9981 ). Note that we have neglected the contribution associated with 
the reduction of the condensation energy cost if a vortex and a flux tube overlap. This leads 
to an energy cost per interesection slightly smaller than that estim ated above (of the order 
E p w 0.1 - 1 MeV) jRuderman et al.lll998l : Isiderv fc Alparl bood ). Vortex motion is thus 
impeded by the flux tubes, that provide an effective pinning barrier unless the vortices have 
enough energy to cut through them. The corresponing pinni ng force pe r unit length of a 
vortex has been estimated to be f p ~ 3 x lO 15 !?^ 2 dyn cm -1 ( Link 2003 ). and is balanced 
by the Magnus force for a critical relative velocity of w c ~ 5 x lO 3 !?^ 2 cm s_1 - This leads 
to a critical lag (at a radius of 10 km) Afi c « 5 x lO -3 !?^ 2 ra d, where we have assumed an 
average density for the core of p = 3 x 10 14 g cm -3 . Given the large value of the critical 
lag, comparable to what could be built up in-between Vela glitches, a substantial part of the 
vorticity in the core could be pinned. 

To account for this effect we shall assume that a fraction of the vorticity in the core is 
pinned and does not contribute to the angular momentum stored in the vortex sheet. This is 
equivalent to assuming that all the vorticity within a radius R v = rjRb is frozen, with Rb the 
radius of the base of the crust and r\ a free parameter. We thus define a fraction of pinned 
vorticity in the core as £ = (-R^)/(^ max )- The total number of vortices in the vortex sheet 
before the glitch scales accordingly: 



;i-o 



2vr 



/v 



(2) 



By using (T5]) in (CD) we can obtain the angular momentum exchanged during the glitch and 
by fitting the size of a glitch, AQ g \, we can derive the coupled fraction of superfluid 



Y, 



1 

Q(i-0 



" AL 



+ Q-1 



(3) 



tot 
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The instantaneous step in the frequency derivative then follows from 

Aflgl Q(l-Q(l-Y gl ) 

^oo i - g[i - (i - era 



Results 



In order to compare our results with observations we consider the case of the Vela 
pulsar. The Vela (PSR B0833-45 or PSR J0835-4510) has a spin frequency v » 11.19 Hz 
and spin-down rate z> w —1.55 x 10~ n Hz s -1 . Giant glitches are observed roughly every 
thousand days and have relative frequency jumps of the order AQ/Q 10~ 6 . The spin-up is 
instantaneous to the a ccuracy of the data, with upper limi ts of 40s for the rise time obtained 
from the 2000 glitch JPodson. McCulloch fc Lewisl l2002h and of 30 s for the 2004 glitch, 
although this limit was less significant ( jDodson et all 120071 )). The glitch is usually fitted to 
a model consisting of permanent steps in the frequency and frequency derivative and a series 
of transient terms. It is well known that to fit th e data at least three are required, with 
decay timescales that range from months to hours (|Flanaganlll996l ). Recent observations of 
the 2000 and 2004 glitch have shown that an additional term is required on short timescales, 
with a decay time of approximately a minute. Given that the Vela 2000 glitch provides 
the most robust observational results, we shall compare the expression in (j3j) to the step 
in frequency derivative associated with the short timescale (1 minute) after the Vela 2000 
glitch, which we assume is a reasonable approximation to the instantaneous post-glitch step 
in the spin down rate. The para meter Y g i is obtained from (|3l) by fi tting to the Vela 2000 
glitch size of AVt/Vt = 2.2 x 10" 6 JDodson. McCulloch fc Lewisl l2002h . 



In figure (pQ) we show the results for varying values of £, for both SLy and GM1. The 
parameter (3 encodes the reduction of the pairing gap due to polarization effects in the 
neutron medium. Recent cal culations suggest tha t polarization reduces the gap and shifts the 
maximum to lower densities (jGandolfi et ajj|2009j ), an effect which in our setting corresponds 
to the value (3 ~ 3, while (3 = 1 corresponds to a bare particle approximation. The horizontal 
lines show the region that is allowed by the measurements of the step in frequency derivative 
of the Vela 2000 glitch. It is obvious that most equations of state and proton fractions can 
mat ch this value if all vo rticity in the core is free, as was also found by lSeveso et all (l2012at ) 
and lHaskell et all (12012al ). although we note that for the more realistic case of (3 = 3 and three 
body cases included in the calculation of x p , the stiffer equation of state is clearly favored. 
We now compare this to the case in which part of the core is in a type II superconducting 
state and part of the vorticity is pinned. As we can see from figure (PQ), as the parameter £ 
increases it becomes increasingly harder to fit the observed values of AQ and in most cases 
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M[M S ] M[M S ] 



Fig. 1. — We plot the value of the fractional step in frequency derivative Ati/Cl for varying 
values of the fraction of pinned vorticity in the core, £ for SLy and GM1. The parameter (3 
encodes in - mediu m polarization effects, as described in the text. For x p we use the results of 



Zuo et al.l (120041 ). both those obtained with two-body interactions (case a) and with three- 



body forces (case b). The end of the curves in these two cases corresponds to the point after 
which we can no longer find a reasonable physical solution. The horizontal line represents 
the measured value for the Vela 2000 glitch, and the thin lines are respectively la and 2a 
deviations. It is clear that in general both EOSs and models for the proton fraction are 
compatible with free vorticity in the core. As we increase the pinned fraction however it 
becomes increasingly more difficult to fit the data, and for £ > 0.5 no solution can be found 
at the la level. 
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this is only possible for a very restricted interval of masses. In general one cannot obtain a 
physically reasonable fit if more than half of the vorticity in the core is pinned. This points 
to the conclusion that the vortex/flux tube interaction is weaker than previously assumed 
and th at most of the vorticity in the core is, in fact, free. This conclusion is compatible with 
that of lHaskell et al.l ( ]2012al ). who found that a weak coupling between the superfluid and 
normal component in the core (as would be the case if most of the vortices in the core are 
pinned) does not allow to fit the shorter post- glitch relaxation time- scales of the Vela. Note 
that the conclusions of this Letter and those of lHaskell et al.l (j2012al ) are derived in different 
methods (in t his case calcu l ating t he exchange of angular momentum in a static model, 
in the case of lHaskell et al.l ( )2012al ) by fitting the post-glitch relaxation with a dynamical 
m ultifluid model) and are thus independent, sa ve for the use of the pinning forces calculated 



in 



Grill et al.l (120121 1: iGrill fc Pizzocherol f l2012l ). 



Conclusions 



In this Letter we have extended the "snowplow" model of iPizzocherol (120111 ) to account 
for the possibility that part of the vorticity in the core may be pinned due to the interaction 
between vortices and flux tubes. We fit the step in frequency and in frequency derivative of 
the Vela 2000 glitch to obtain constraints on the pinned fraction of vortices in the core and 
in general find that both quantities cannot be fitted for reasonable physical parameters if the 
pinned fraction is larger than 50%. Although we do not deal with the microphysical details 
of the vortex dynamics in the core, our conclusions are quite general. The only quantity 
that is needed to evaluate the angular momentum that is exchanged during a glitch is, in 
fact, the number of vortices that are stored close to the peak of the pinning potential in the 
crust. As long as the excess vorticity of the core can be transferred to the equatorial strong 
pinning region in-between glitches the details of the vortex motion are not influential. 

The general conclusion is that either most of the core is in a type I superconducting 
state (and the vortex pinning is negligible (jSedrakianl 120051 ) ) . or that the vor tex/flux tube 



i nterac tion is weaker than previously thought. This is the same conclusion that lHaskell et al. 



(I2012al ) come to after fitting the post-glitch short-term relaxation of Vela glitches with a 
hydro dynamical model. If such an conclusion is confirmed it would have serious implications 



20 12t iHo et al.ll201ll ; lHaskell et al. 



also f or NS precession (iLinki 12003 ) and f or G W emission (lHaskell et al.l 120081 ; lLander et al. 



2012b|). Note that on a microphysical level it is very likely 



that the interaction between vortices and flux tubes is weaker than the estimates presented 
here. These estimates are upper limits on the strength of the pinning force, as they do not 
account for the finite rigidity of vortices, which could lead to a reduction of a factor 100-1000 
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flLinkl l20 



2b; Grill & 



izzocheroll2012l ; iGrill et al.ll2012[ ). as rece nt theoretical estimates also 



confirm (ISeveso et al.ll2012bf ). Furthermore recent calculations (lBabaevll2009l ) suggest that 
in the presence of strong entrainment or gapped £~ hyperons in the crust the interaction 
between flux tubes and vortices will be significantly weaker, and even in the presence of 
pinni ng the superfluid may be coupled to the crust on short timescales ( iSidery fc Alpar 
20091 ). It should also be pointed out that if flux tubes are able to move out with the neutron 
vortices on the interglit ch timescale, this could le ad to a potential barrier at the crust core 
int erface and a glitch (ISedrakian fc Corded Il999l ). However recent quantitative estimates 
by iGlampedakis fc Anderssonl (120111 ) indicate that in the pinning regime vortices and flux 
tubes can be considered essentially immobile on the inter-glitch timescale, as we assume in 
the present model. 

Finally note that we have assumed straight vortices that cross the core. Superfluid 
turbulence is, however, a wel 
to also occur in NS s (see e.g. 



known phenomenon in labo r atory superfluids and is expected 



Andersson. Sidery fc Comerl (120071 )). especially in the presence 



of strong pinning (jLinki l2012bl Jal). In this case the vortices will form a turbulent tangle, 
leading to a weaker inter action between the components and longer c oupling; timescales in 
the outer core and crust (IPeralta et al.l 120061 ; iPeralta fc Melatod 120091 ) , which will however 
couple the components on long inter-glitch timescales. Vortex pinning in the core would 
then be a transitory phenomenon and most of the vorticity would be free on long timescales 
(see also Glampedakis &: Andersson (j2011 )) . The transition between a laminar and turbu- 



lent flow could however have a strong impact on the glitch mechanism (IPeralta et al.l 12006 



Melatos fc Peraltal 120071 : IPeralta &: Melatod |2009| ) and the definition of a pinning force per 
unit length is significantly more complex in a vortex tangle. We plan to tackle this issue in 
future work. 
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